March 1994 



EQUATIONS OF GRAVITATIONAL 
INSTABILITY ARE NON-LOCAL 



Lev Kofman 

Institute for Astronomy, University of Hawaii, Honolulu, HI 96822 

and 

Dmitry Pogosyan 

CITA, University of Toronto, Toronto, ON M5S lA7,Canada 



Submitted to: The Astrophysical Journal 



1 



ABSTRACT 



Few recent generations of cosmologists have solved non-local newtonian equations of 
the gravitational instability in an expanding universe. In this approach pancaking is the 
predominant form of first collapsing objects. Relativistic counterparts of these equations 
contain the electric and magnetic parts of the Weyl tensor. In the linear theory the 
magnetic part is associated with gravitational waves. If the magnetic part is ignored, 
then the newtonian limit of the relativistic equations is reduced to the closed set of the 
local Lagrangian equations. Recently this fact drew much attention since the gravitational 
instability in that form would greatly simplify the study of cosmic structure formation. In 
particular, the filamentary structure of collapsing is predicted. In this paper we resolve the 
contradiction between the newtonian theory and relativistic version adopted in some recent 
papers. We show that dropping the magnetic part from the basic relativistic equations is 
incorrect. The correct newtonian limit is derived by the 1 / c-expansion of the GR equations 
and the Bianchi identities for the Weyl tensor. The last ones begin with ~ 1/c 3 order, 
therefore one must take into account the magnetic part in the post newtonian order ~ 1/c 3 , 
which contains non-local terms, related to the non-local gravitational interaction. For 
the first time we rigorously show that the basic GR equations with the magnetic part 
are reduced precisely to the canonic newtonian non-local equations. Thus, the correct 
treatment of the relativistic version of the gravitational instability resurrects the canonic 
picture of the structure formation. 

Subject headings: cosmology: theory — large-scale structure of the universe 
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1. INTRODUCTION 



The cosmological gravitational potential related to the large scale structure within the 
horizon is small in the sense of General Relativity, (^<c 2 , where c is the speed of light. 
Motion of matter is also relatively slow, « «c. Hence within the cosmological horizon 
one can use the Newtonian gravity to describe the gravitational instability of matter in 
an expanding universe (Zel'dovich and Novikov 1983, Peebles 1980). Newtonian theory 
is also applicable in the most of cases of gravitationally maintained astronomical systems 
( Binney and Tremaine 1987). Newtonian theory formally can be derived from the 1/c 
decomposition of the Einstein equations in the limit c — > oo. The Poisson equation, linear 
with respect to 0, arises in 1/c 2 order of the GR equations for gravitational field. 

The set of equations describing the self-consistent evolution of the pressureless matter 
of density p with newtonian gravity is non-linear in terms of p, v and <fi. It is clear, 
intuitively, that the instability of the media due to the gravity in the newtonian theory 
is also non-local, since the gravitational potential at a given point is a superposition of 
contribution over all distributed masses, as textbooks teach us (Zel'dovich and Novikov 
1983, Peebles 1980, Tremaine and Binney 1985). The dynamics of the non-linear non-local 
gravitational instability is very complicated even in the reliable newtonian theory, and 
requires the N-body simulations of the non-local gravitational interactions and motions 
of matter. Therefore any analytical advantage would be greatly appreciated. There are 
very few successful analytical approximations to the solution of the basic equations. One 
of the best known in the cosmological context is the generic Zel'dovich approximation 
(ZA, Zel'dovich 1970). ZA describes the displacement of particles with the initial velocity 
in the quasi-linear regime, which is realized at early stages of non-linear evolution, or at 
sufficiently large scales, where evolution is also quasilinear. ZA is local, in a sense that 
particle evolution upon gravitational forces truncated along ZA, depends on its Lagrangian 
position only. ZA is an approximation and fails after the particle orbits' crossing. ZA 
quantitatively predicts that the first non-linear objects in any cosmological scenario are 
two-dimensional pancakes, which then evolve into a cellular structure with the most mass 
in the clumps or filaments, depending on the spectra. Further, non-linear evolution looks 
like complicated superposition of hierarchical clustering and pancaking over a vast range 
of scales, with qualitative characteristics crucially depending on the model (see, e.g., 
Shandarin & Zel'dovich 1989, Kofman et al 1992 for theoretical outlook, Cen & Ostriker 
1992, Bertschinger & Gelb 1991 for numerical outlook). 

One would expect that relativistic generalization of the self- consistent matter evolution 
in an expanding universe is much more complicated because of the additional non-linearities 
in the Einstein equations. In 1971 Ellis reviewed an elegant form of GR equations together 
with the Bianchi identities of the self-consistent matter evolution, suggested earlier by 
Triimper 1967. In GR, one can split the Riemann tensor, Rikim, which entirely describes 
the structure of gravitational field, into the Ricci tensor defined by the Einstein 
equations, and the Weyl tensor Ciki m (e.g., Landau & Lifshitz 1980). One can rewrite the 
Bianchi identities in the form rather reminding of the Maxwell equations (e.g., Hawking 
and Ellis 1973). Additionally, one can split the Weyl tensor into the "electric" part Eij, and 
the "magnetic" part Hij, due to some of their similarity to electrodynamical counterparts. 
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In general relativity 00 component of the Einstein Eqs. - relativistic counterpart of 
the Poisson equation - plays the role of the constraint equation. The constraint equation 
is automatically resolved for an arbitrary time moment t > to if it is resolved once at the 
initial hypersurface t = to, and the evolution equations are satisfied. Note that in the 
Newtonian limit the Poisson equation has to be resolved at each time step t. Ellis gave a 
set of evolution equations for E^ , and Hij . Equation for the Lagrangian derivative of E^ 
contain their algebraic combination, except one non-local term - the space derivatives of 
the magnetic part Hij-k- 

For the electric part in the c — > oo limit we put E^ = -^E^ + -^e\J N ^ + O(^). 
The Newtonian limit of the electric part is reduced to the tidal forces 

=V Xa V x ,<j ) -U a pVl<j ) , (1) 
e[J NS> is the first post-newtonian term. Respectively, for the magnetic part 

*«=?4* + <>(?). (2) 

The magnetic part in the linear theory contains contributions from the vector and 
tensor modes only. Vector modes are irrelevant for cosmological gravitational instability, 
contribution from cosmological gravitational waves is negligible, they are generated by 
matter in 1/c 5 order (in the Lagrangian) only (e.g. Landau & Lifshitz 1980). Therefore 
we can put the newtonian term = 0, and H^ N ^ ^ is the first non-vanishing 

post-newtonian term. There are controversies in the literature at this point. For instance, 
it is incorrect to say that the magnetic part is non-zero in the newtonian limit in the absence 
of vector and tensor modes. It is non-zero in the first post-newtonian approximation only. 

Since the magnetic part is vanishing in the Newtonian limit, there is a temptation 
to put it always equal to zero in the evolution equations (Barnes & Rowlingson 1989, 
Mataresse et all 1993). Then one could obtain the closed set of the Lagrangian local (but 
non-linear) equations for Eij,v and p. It looks surprisingly simpler than the non-local 
equations of the Newtonian theory, since one can deal with the set of ordinary differential 
equations. 

Recently this fact drew much attention in connection with possible application to 
the dynamics of gravitational instability beyond the linear theory (Mataresse et all 1993, 
Croudace et al 1994, Bertschinger & Jain 1994). Some far-reaching conclusions were 
made based on this interpretation of the relativistic equations. For the closed set of the 
Lagrangian equations which essentially are the ordinary differential equations, the phase 
diagram can be drawn. The instability of collapse in this theory was found (Mataresse 
et all 1993). Bertschinger & Jain 1994 concluded that typical collapsed configuration are 
strongly prolated (one dimensional) filaments rather than Zel'dovich pancakes. 

The implication of these new equations for the Newtonian theory remained unclear. 
There is a clear contradiction with the canonic newtonian picture of the gravitational 
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instability in cosmology. Bertschinger & Jain 1994 cautiously noted the lack of any 
Newtonian derivation of the set of the closed Lagrangian equations in this approach. 
However, the confrontation of views has not been resolved. In this paper we give the 
resolution of this contradiction. In Sec 2 we recall the basic equations in the Newtonian 
theory, and prove that they are non-local for physically interested values, like the tidal 
forces. We derive the newtonian equations for the tidal force (1). It will be a convenient 
form in order to compare them later with the relativistic calculations. In Sec 3 we rigorously 
implement the transition to the newtonian limit of the relativistic equations reviewed by 
Ellis 1971. We show that the equation for En- begins with 1/c 3 order, therefore the term 
with the magnetic part Hij-k involving in this equation, has to be calculated in 1/c 3 order, 
i.e. beyond the linear newtonian limit. To calculate the components of the Riemann tensor 
up to 1/ c 3 terms, we have to take into account the first post- newtonian corrections in the 
Einstein equations - the classical problem, considered many years ago by Einstein, Infeld 
and Hoffman (1938). Using this result, we derive the equations reviewed by Ellis 1971 in 
the first non- vanishing order of 1/c decomposition, which turns out to be of /c 3 order. 
We demonstrate, that the resulting equations are precisely reduced to the results of the 
newtonian theory. In the inverse order, from the newtonian equations we can rederive the 
GR version of the basic equations in the first non- vanishing order of 1/c decomposition, 
and show that the term related to the magnetic part inevitably arises. Thus, dropping the 
magnetic part from the basic GR equations would have a dramatic effect similar to the 
dropping of the magnetic part from the Maxwell equations. 

Among the published works so far, our view of the problem is closest to that of 
Matarrese et al. 1994, who also claimed that putting = cannot be exact. They 

used the second order in the cosmological perturbations series, and found gravitational 
waves beyond of the horizon, whose nature is not clear to us. We use 1/c-series up to 1/c 3 
order in metric gik, where gravitational waves are not yet generated and do not interact 
with matter, and therefore need not be considered. It would be interesting to perform 
both decompositions at once. 

We conclude that the correct basic equations of gravitational instability are 
intrinsically non-local. In Sec. 4 we additionally discuss the Zel'dovich solution, which 
is the local approximations to the cosmological gravitational instability. We give the 
explicit solution for the tidal forces in the ZA. In the conclusive Section 5 we argue that 
two-dimensional pancaking is the generic form of collapsing in the single stream regime. 



2. NEWTONIAN TIDAL FORCES ARE NON-LOCAL 

2.1 Newtonian equations via potentials 

In the epoch of the LSS formation, most of the mass is in the form of dark matter 
of relic origin (for instance, the Cold Dark Matter), with no pressure; nonlinear Large 
Scale Structure is originated by gravitational instability of small initial fluctuations. Let 
x, v = a^j|, p(t, x) and <f>(t, x) be, respectively, the comoving coordinates, peculiar velocity 
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and density of dark matter (neglecting baryons), and peculiar Newtonian gravitational 
potential. In the Newtonian theory, the motion of dark matter before the particle orbits 
cross obeys a nonlinear system of equations 

^+3^p+-V x (pv) = 0, (3) 
ot a a 

dv 1 _ a _ 1 _ , , t . 

— + -(v-V x )v + -v = —V x <j), 4 
ot a a a 

V x 2 c/> = 47rGa 2 (p-p), (5) 

where a(t) is a scalar factor of an expanding Universe and p is a mean density. These 
equations are valid in the single stream regime, and admit the evident generalization in 
the regions of multiple streams. 

For the sake of simplicity, we assume that inflation produces a flat universe = 1, 
and A = 0. The growing mode of adiabatic perturbations D(t) in the Einstein-de Sitter 
universe is D(t) = a(t) oc t 2 / 3 . However, all results reported here are valid for an arbitrary 
background cosmology. The solution of the Poisson equation (5) is non-local, therefore the 
motion of the fluid element upon the gravitational forces in the right hand side of the Euler 
equation (4) is also non-local, i.e is not entirely defined by the masses along the particle 
trajectory. 

One can also conclude that the derivatives of the gravitational force - the tidal force 
- given, for instance, by equation (1), are also non-local. However it would be useful to 
write down non-local equations for the newtonian tidal forces, or some other combinations 
of its further derivatives (Kofman 1993), especially since there are opposite statements in 
the literature (e.g. Lachieze-rey 1993). For this purpose let us perform further analysis 
of the basic equations (3)-(5). It is convenient to use the growing solution D(t) = a(t) 
as a new time variable instead of t (Gurbatov et al 1989, Kofman 1991), and introduce a 
comoving velocity u = = v/aa in respect with this time variable. Then eq. (4) acquires 
the form 

+ (u- V x )u = (ad)" 2 ( (3a 2 a)u - V x </> ) (6) 



da 

Let $ be velocity potential so u = V x $. We introduce the combination A = (|i7 2 a 3 ) -1 = 
— (3aa 2 ) -1 = (47rGpo) _1 5 which does not depend on time in the matter-dominated 
Einstein-de Sitter Universe. Then from eq. (6) we find a general relation between the 
velocity potential and the Newtonian gravitational potential (Kofman and Shandarin 
1988): 

^ + k(V**) 2 = -£(* + VI (7) 



da 2 V ' 2a 

This equation is often referred to as the Bernoully equation. This form is very convenient 
to discuss the Zel'dovich approximation, in which the right hand side is vanishing (Kofman 
1991). Substituting density p from the right hand side of eq. (5) into eq. (3), we find the 
second equation linking scalar potentials <j> and $ 

d 

-{aV x c4>) + (A' 1 + aV 2 ^)V x «$ 



da 



0. (8) 
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Thus we obtain two equations (7), (8) for two scalar fields $ and 4> m the single stream 
regime. Note, that from (7) one can express <fi via $ and its derivatives, and substitute it 
in this form into (8). Then one gets the single non-linear partial differential equation for 
$ only. 



2.2. Tidal forces in Newtonian theory 



The structure of equation (8) is V x • £ = 0, where the vector is the expression 
in the square brackets in eq. (8). Using the Helmholtz theorem of the vector calculus, we 
find 



£(£) = V x x 



47T / \X — X 1 



(9) 



Substituting the expression for from the square bracket of eq. (8), in eq. (9), after 
some vector algebra, we end up with the formula 



(10) 



In the Newtonian theory we can construct the traceless symmetric tensors 
along with the traceless tidal forces tensor (1). Taking derivatives V X «V^ from eq. (7), 



we can derive the equation for the Lagrangian time derivative of along the trajectory 

- _ A (J N ) + AE (N) ] (11) 



Taking derivative V x /3 of eq. (10), we derive the equation for the evolution of the tidal 
forces along of trajectory 



d 3 x 



a/3 



47T 



\x — X 



-(V X ;$)(V^) 



(12) 



At this point we rewrite eq. (12) back in the coordinate time and physical variables 
to compare it with the result we will derive in Sec. 3 from the relativistic equations: 



d 



(N) a ( N )_ 



4nG 



dt af3 



aV x (pv)d a/3 = aGV x pV x o 



J 



d 3 x' „ , . 
- — ^x'Kpv) 



\x — X' 



(13) 



7 



Note the non-local integral term selected to the right hand side. This non-local term 
cannot be reduced to the combination of the local terms. It prevents us from obtaining a 
closed set of Lagrangian equations in the rigorous Newtonian theory. 



3. TIDAL FORCES IN GENERAL RELATIVITY ARE NON-LOCAL 

3.1. Notations and GR equations 

In General Relativity the gravitational field, described by the Riemann tensor Rikim-, 
is determined by the Einstein equations 

= (V - \6}T,y) (14) 

and the Bianchi identities 

2C lklm . m = R H > k - R lk ' 1 + ^ (g lk R> 1 - g li R' k ) (15) 

Here R k = R m% mk is the Ricci tensor, R = R m m is a curvature, Cikim is the Weyl tensor 
and () ; is a covariant derivative. Latin letters run from to 3 and Greek letters run from 
1 to 3. The signature of metrics we use is ( — h ++)• For a pressureless matter treated 
as the perfect fluid, energy- momentum tensor is given by T tk = eu l u k : e is the energy 
density and u % is four- velocity of the fluid, u l Ui = — 1. Equations of motion of the matter 
contained in the Einstein equations are 

T k k , % = (16) 

One can further decompose the Weyl tensor into the "electric" part E ik = Cnk m u l u m 
and the magnetic = \r\u mv C mvkq u l u q (where rjikim is fully antisymmetric four-tensor): 

Cikim = (VikpqVlmrs "H 9ikpq9lmrs)' l ^U (j]ikpq9lmrs ~t~ Qikpq^lmr s)^ ^ ? 

where 

Qiklm = 9il9km - Qimdkl- (17) 

Then the Bianchi identities lead to the "Maxwell-like" equations on E ik and H ik (Triimper 
1967). From them of the most interest is the evolution equation on E^ which we write 
down for the case of pressureless irrotational fluid 

p l n p k l u rn E nl . m +p lk a lm E lm + GE ik - 3E l{l a k) l + Pi {l i] k)r sm u r H l s . m = -^-ea lk . (18) 

Here p lk = g lk + u l u k is the projection tensor. The velocity derivative tensor is presented 
as the sum of its trace 0, symmetric traceless shear a ik and antisymmetric vorticity Ui k 
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as PuPk m u l -m = Cik + \®hik + u>ik- In our case lo^ = 0. Parentheses on the pair of 
indexes denote symmetrization, while brackets later in the text - antisymmetrization. The 
first term in (18) is the Lagrangian derivative of the electric part along of the trajectory, 
other terms are algebraic combinations of the local tensors, except one non-local term, 
containing the space derivatives of the magnetic part H l s . m . The rest of the equations 
that follow from the Bianchi identity the reader can find in Ellis 1971. Note that there are 
two local equations for and © in the absence of vector and tensor modes. 

3.2. | expansion of the Weyl tensor for the gravitational instability in an expanding 
universe 

Now we shall analyse the newtonian limit for the GR equations that describe the 
motion of the self-gravitating pressureless irrotational fluid in an expanding Universe and, 
specifically, the role of eq.(18) for the electric part of the Weyl tensor in this limit. 

To describe the gravitational instability on the uniform expanding cosmological 
background, let us write the metrics as 

ds 2 = -a 2 (r) [c 2 (l + h 00 )dT 2 - 2h 0ce cdrdx a - (S a/3 - h af3 )dx a dx f3 ] , (19) 

t is the conformal time, for the sake of simplicity considering the flat space geometry. 

Newtonian limit of the GR equations formally can be obtained by taking the limit 
1/c — > 0. We will need the decomposition of the metrics hij up to first post-newtonian 
correction ~ 1/c 3 . The quantities describing matter behave themselves as e — > pc 2 , 
u° — > a -1 , uq — > —a and au a = c~ 1 v a , where p is the mass density and v a - physical 
velocity. Hence the right-hand side of eq.(14) is of order of 1/c 2 and, as we are interested 
in self-consistent gravitational field produced by matter, the expansion of hik also starts 
from 1/c 2 terms. Substitution of the metric in form (19) shows that when 1/c — > only 
/loo remains present in the equations of motion (16): 

= (20) 
= -\a 2 h^ a (21) 

Here the prime denotes derivative with respect to conformal time r while ordinary spatial 
derivatives are denoted by () ja . The correspondence with standard 3D notation we used 
in previous Section is: f = i£ = £()', V x « = ^ = a 2 ^- = (), Q . Together with (00) 
component of the Einstein equations (14) which determine hoo: 

\a 2 h m ^= A -^-a 2 {p-p) (22) 

the equations of motion, thus, form a closed set of equations for /i o 5 P and v a , identical 
to the Newtonian ones (3)-(5) if one identify hoo with gravitational potential hoo = 2^/c 2 . 
Next post-newtonian correction in hoo is l/c 4 -order. 



a 



a 



9 



The remaining components of the Einstein equations simply determine ho a and h a p 
that are consistent with the Newtonian limit through the known quantities 0, p, v a . Using 
well-known formulas for post-newtonian metric decomposition, e.g. from Landau & Lifshitz 
(1980, §106), and generalizing them including the expansion of the universe, we have 



Kp = h 00 5 a p + O [ \) = \2(/)5 a (3 + i^j 



1 1 

C 3 4"7T 



/^1( 167rG ^ + ^ (a4 ^ }/ )- 



(23) 
(24) 



In the last expression the auxiliary condition a 2 ho a ,ce = \(a 2 h a a ) is assumed. Let us 
note that the quantity ho a is already non-local. 

The expressions (23)-(24) are sufficient to write down the Weyl tensor Cikim up to 
the order 1/c 3 and to determine the first non- vanishing terms for E ik and H ik : 



Co a o(3 — a 2 E a [3 
Coaf3j = a 2 



^05, [7 ) gp]a 



Cap-ys = E aj gp$ + Ep S g ai - E aS gpj - Ep^g aS 



1 



, ,1) _ 1 UN) 



H a p = VaS^Cyopou 6 + -Vc e o' yS C 1 spT 



,u 



(*) 


(25) 




(26) 


(*) 


(27) 


uo 


(28) 


(*)■ 


(29) 



The remaining (00) and (0a) components of E ik and H ik tensors are of a higher order in 
1/c. 

Although Newtonian equations (20)-(22) arise as ~ 1/c 2 terms in the expansion of 
relativistic equations, (and to obtain them one does not need to consider any higher order 
quantities), it is no surprise, however, that one can construct such combinations of GR 
equations which do not contain 1/c 2 terms at all. In this case, to reproduce correct limit 
1/c — > of such combinations it is necessary to carry out the calculations beyond 1/c 2 
order. At the first glance one could think that the new relations between Newtonian 
quantities could arise as a result, but, of course, it is evident that all these new equations 
must appear to be some combination of basic set (20)-(22) in order for GR to be compatible 
with the Newtonian theory. The evolution equation (18) for E^ just provides us one such 
example. Let us fix free indexes in eq.(18) to be spatial ones (a/3). Using expressions 
(28)-(29) and the Newtonian limits for the shear aa a p = \cr^p and the expansion a© = 
-@( N \ we get in leading order 



1 



J ap 



E^V'+-E^+E^ v-r 
a 



J a(3 



(n) (N)-y 

7 (a a /3) 



+ AnGpa 



(TV) 
a/3 



-V(po sm H a)s , m (30) 
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The left hand side of this equation is of order 1/c 3 . The magnetic part of the Weyl tensor 
at the right side is non- vanishing in this order, according to (29). Hence, it does not 
drop out from eq.(30). Therefore we conclude that despite the fact that is of higher 
order of smallness than it cannot be set to zero in eq.(18) on the basis of considering 
Newtonian limit. As we show later evolution equation for E^ derived from GR remains 
to be nonlocal. 

In the next section we proceed to prove that Bianchi identities and, consequently, 
eqs.(18),(30) are completely consistent with standard Newtonian equations (3-5). 



3.3 Reduction of Bianchi identities into Newtonian equations 



The straightforward way to investigate eq.(30) in the Newtonian limit is to insert 
expressions (28), (29) for E^ and into (30). The following algebraic computations are 
rather lengthy and given in the Appendix. It is more convenient to deal directly with the 
Bianchi identities utilizing the components of the Weyl tensor given by (25)-(27). To get 
rid of Rik terms we rewrite eq.(15) using the Einstein equations as 



c-» 



(31) 



The principal combinations of indexes (ikl) are (OctO), (a0(3) and (a(3^). In the leading 
order in 1/c, combination (OaO) straightforwardly corresponds to the derivative of the 
Poisson equation 



E, 



a(3 



8tvG 
"3c2" 



p >o (^ >7 -47rGp) |O = 



(32) 



The (af3"f) combination leads to the equation for the time evolution of the magnetic part 
Hi^. However, in the Newtonian limit all terms, containing H^, appear to be of higher 
order than the rest of the equation which gives us again 



(33) 



The fact that the relativistic evolutionary equation for reduces in the Newtonian limit 
to an identity that does not include time derivative of suggests that magnetic part 
of the Weyl tensor does not have an independent dynamical meaning in the Newtonian 
theory. 

The last combination (oeQ/3) gives rise, in particular, to the evolution equation for E^ 
and is slightly more complicated. Substituting (26) into (31), we get 



a' 
+ -E, 

a 



a(3 



a i S 



4nG 



,/3]o 



(34) 
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Noticing that the antisymmetric term in the squared brackets the left hand side is [ ...] = 



^ap + 7T0,a/3 - we finally have 

which is essentially 



d_ o!_ 

dr a 



(A(f>-47rGa 2 (p- p)) = 0. 



(35) 



(36) 



If we return to eq. (34) and substitute there ho a from eq. (24), we arrive at the form 



8 (AH a' (AT) 4tvG 2 



dr ^ 



\x — x 



T\P V 



(37) 



,7a/3 



which is completely identical to eq. (13) derived from the Newtonian theory, but only 
written in the relativistic notations using r and g a p. Thus, we prove that the relativistic 
equations for the gravitational instability of the scalar perturbations in an expanded 
universe are reduced to the newtonian equations, which are non-local. 



4. LOCAL APPROXIMATIONS 

As we demonstrated in the previous sections, rigorous Newtonian equations of the 
gravitational instability in an expanded universe are intrinsically non-local. However, this 
does not exclude some useful intermediate approximations which might be local. Famous 
Zel'dovich approximation (Zel'dovich 1970) is an example of the approximated description 
of the gravitational instability, which turns out to be a local one. 

It is convenient to use the equation of motion in the Lagrangian description. In such a 
case the dynamics is described by the displacement, S(q, t), of each particle from its initial 
position q. Its current Eulerian comoving position, x, is then given by 

x = q + S(q,t). (38) 

For ID case the mildly non-linear dynamics is quite simple. The reason is that the 
force exerted by a density perturbation over a given particle is independent of its distance 
to the particle. Therefore, before any shell crossing, the displacement of each particle 
depends on its Lagrangian position q only. Then the displacement field for the growing 
mode can be factorized S(q,t) = W^(q) D(t). In ID case the ZA is then identical to the 
exact dynamics. In relation (2) the displacement factorizes in a spatial function, V^(q), 
which depends on the initial conditions, and a universal time dependent growing mode 
D(t). Using equations (7), (8) for the potentials 4> and $, one can rigorously derive in ID 
case that 

(f ) = -A- 1 §, (39) 
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and then the dynamical equation (7) for the velocity potential $ is reduced to the 
"shortened" equation without the right hand side: 



f + i(V»*) 2 = 0. (40) 

The 3D generalization of the form (38) of the displacement is the Zel'dovich 
approximation (ZA). The 3D generalization of the equations (39) and (40) is the Zel'dovich 
approximation expressed via velocity potential Note the simple solution of equation 
(40): 

d>(a,f) = MO) + = *o(9) + \A (41) 

The solution (47) describes the deformations of 3D-hypersurface of the potential and 
is valid until the formation of folds of $-hypersurface, which corresponds to caustics 
(pancakes) . 

In this section we consider the equations for the tidal forces in ZA, show that they 
are local, and find its solution. Let us introduce the tensor of the velocity derivatives 
W a p(x,t) = —VxocVp in the Eulerian space. For the potential motion it is reduced to 

W aP = -V xa V X/3 $. Note also that W a p = -a^J - \8 a p& z \ Here = V*$. Let A< 
be its eigenvalues. The field of the W a p(t)-tensoT evolves in time, its initial value (in the 
Lagrangian space) coincides with the Lagrangian deformation tensor D a p = — V q <* V q p<&o, 
Aoa are its eigenvalues. We will call the S^-tensor the Eulerian deformation tensor. The 
tensor of tidal forces is the ZA, with help of is equation (39), is 

E^p = ~j(w a p- 1 -5 a pW^, (42) 

Taking space derivatives of the equation (40) and using its definition, we derive the equation 
for the tidal forces in the ZA. 

IL e (Z) + — e( z ) E (Z) + E^ T a {Z) - —5 ■ B( z )/r( z h T - (43) 

This equation is simpler than equation (12). Actually, in ID case equation (12) reduced 
to equation (43). 

It is easy to construct the solution of equation (43) for the tidal forces in the ZA Using 
the definition of W a p and its eigenvalues X a (t), one can instantly find from (41) how they 
are related to their initial values Ao a : 

■ (irki- (44) 

Then the principal radii of curvature R a = A" 1 linearly increase with a(t): R a {t) = 
Roa + d(t). It means that the $-hypersurface obeys the rule of the evolution, which is 
similar to the Huygens principle of geometric optics for the wave front. 
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Substituting (44) into (42), we conclude, that the principal axis of the tidal forces 
tensor are not rotating in the single stream ZA, and its eigenvalues ^ are evolving with 
time as 

f*«(t) = ~\ (3A a (£) - (Xi(t) + X 2 (t) + A 3 (£)) . (45) 

Thus, the form (45) is the solution of equation (43). In the work in the progress (Bond et 
al 1994) we are comparing the tidal forces in the ZA against the actual tidal forces. 

Thus, the tidal forces in the ZA are closely related to the deformation of the 3D 
hypersurface $(£), described by the simple Huygens principle. In the ZA the evolved fields 
can be entirely described via the initial eigenvalues A a o-s - ZA is a local approximation. 
It may be interesting to consider the solution of the exact single stream equation (12) in 
terms of decomposition around the form (44), which is the exact solution in ID case. 



5. CONCLUSION 

We implement the 1/c decomposition of the basic GR equations, describing 
gravitational instability of the pressureless matter with the scalar initial perturbations in 
an expanded universe. It allows us to derive the newtonian limit of these equations. These 
equations include the Bianchi identities, which begin from 1/c 3 order. Therefore we have 
to take into account the Weyl tensor up to this order. It requires the first post-newtonian 
corrections to the metric g^. Having borrowed corresponding results from the textbook 
(Landau & Lifshitz 1980), we conclude that the relevant post-newtonian corrections are 
contained in the magnetic part of the Weyl tensor, meanwhile the relevant electric part is 
needed in 1/c 2 order only. Using this post-newtonian magnetic part, we then end up with 
the canonic newtonian equations for the gravitational instability in an expanded universe. 
The reader can find the formula for the magnetic part in the Appendix. 

We should warn the reader about the statements in the literature on the magnetic 
part. In the absence of the vector and tensor modes, it is incorrect to say that the magnetic 
part is non-zero in the newtonian limit 1/c 2 . It is correct to say that it is zero in this limit, 
but it is not zero in the first post-newtonian limit 1/c 3 . 

Thus, we can continue to work with the newtonian equations for the gravitational 
instability. The results reported in the papers exploited the set of local Lagrangian 
equations, are irrelevant for the general case of cosmological gravitational instability. 
However this approach is useful for the special class of the solutions of Einstein equations 
constrained by the condition H ik = (Croudace et al, 1993). 

Another controversial question aroused from the recent discussion regards to the first 
collapsing forms. It was known for a long time that the first collapses occur generally not 
around the density peaks, associated with the maxima of the sum p oc Aoi + A02 + A03, 
but rather around the maxima of the leading component of the deformation tensor Aoi ( 
Shandarin & Klypin 1983). For the high density peaks their positions can be very close. 
Initially underdense regions also can undergo collapse, since maximum Aoi can be realized 



14 



in the region where the sum A i + A 2 + A 3 is negative (Zel'dovich & Novikov 1983). 

Another question is the form of the first collapsing objects. In the general case the 
local density can be obtained by the inverse of the Jacobian of the transformation (38) 
between q and x, so that 



\dx/dq\ |(1- 50(1 -5 2 )(l-5 3 )r 

where S a are the eigenvalues of the deformation tensor V^Sp, S(q,t) is the displacement 
vector. The form of the collapsing objects depends on the rate of decreasing of the factors 
in the denominator (45). In this case the leading eigenvalue S\ reaches the unity alone, the 
two dimensional pancake form first. In the ZA, the displacement vector is expressed via 
the initial velocity, and for Gaussian initial fluctuations 92% of the Lagrangian particles 
undergo collapse via pancaking (Zel'dovich & Novikov 1983). Pancaking is also a typical 
configuration for the collapsing homogeneous ellipsoid (e.g. Peebles 1980). Since in the 
first approximation it is a reliable generic model in the vicinity of collapsing region, 
this is a strong argument for the general character of pancaking. Note that pancakes 
are pronounced for the baryonic or neutrino models, but can have quite complicated 
"fractal" -like structure for the model with power-spectra close to n = — 3 (Kofman et 
al 1992). 
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APPENDIX 



In this Appendix we write down the explicit expression for the term containing 
in eq.(18) in Newtonian limit (leading terms in 1/c expansion). Using this result we 
show that eq. (18) and eq.(34) are equivalent. In 1/c expansion the time derivative and 
spatial components of four -velocity u a are the quantities of c _1 order. The spatial part 
of projection tensor is then p a p = g a p. Fixing indexes (ik) in eq.(18) to be space-like, 
(i — > a, k — > (3), we have 

„ „ sm r tjI c— «x> 1 „ s;m 1 slnpls^i „,d 1 \ ,rn 

Pl(aV/3)r u H s;m > ;10m S (^a) - ^*ras(^ [^a)dnpUlU J 



a 

; m 



U U r y + (C m0r ( a U /3 )U r )' - (C ( a p)rUmU r )' ]. 



— 4!5] p 5^5f5m is used where brackets denote 



At the last step the identity rf r rjikim 
full antisymmetrization, while parentheses on a pair of indexes denote symmetrization. 
Keeping the leading terms in this formulas, we get 

and substituting expressions (25)-(27) for Cikim, we get further for the last expresion 



a 



0(a,/3),7 



- /io 7 ,a,/3' 7 ] + [2^(0,1/^) - 2E a pu 1 - g a pE 1& u s + g l{a Ep )s u 5 ] ' 7 . (Al) 



To proceed further we note that the first term here is equal to 



1 

c 3 a 



C CL 



Next, we use the Einstein equations for -Roa, and decompose four- velocity derivative as 
u a,/3 = o"a/3 + ^©fi'a/J (assuming irrotational pressureless fluid). Keeping in mind the 
Newtonian expression (28) for E a p, it is straightforward to obtain finally 



Pl( a V(3)r Sm U r H 



sm„,r zjl c—*oo 
s:m 



1 fl 



a 

■"-0 + ^2^' a,/3 



AnG . 1 . 
— —p(cr a( 3+-®g a (3)- 



1 



3c 2 



.7 5 



a a z ' ' ' / r- ' v ~" 3 

-E a p a v? + 3£' 7 ( a a 7 |S ) - -Ea/S© - 9apE 1 s<y' yS - (A2) 
Let us now rewrite eq.(18) as 

1 fl 

c 

= -p l{a Vf5)r Sm u r H l s , m . (A3) 

We see that the magnetic part H a p from (A2) contains largely the same terms as a left-hand 
side of equation (A3), which therefore cancel out. What is left takes the form (6 



\- E a(3 + ~2 E ap ) + E^p^u 1 + g a pa lS E j5 + QE af3 - 3E 7(a a /3 ) 7 + ^^-pa a p 

V CL CL J C 



u 



1 



3c 3 



(9c 



, 7 ' 7 ) / + — ^a/30, 7 ' 7 

a 



1 

= — ~a 
3c 2 



This expression clearly can be reduced to the derivatives of the of Poisson equation: 
1 



;9a(3 



d_ a' vl d 
dr a dx 1 



\S7<j)-AiiGa 2 {p- p)] = 0. 



(AA) 



{AS) 
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